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In this paper, we first define the notion of viscosity solution for tire following system of partial 
differential equations involving a subdifferential operator: 
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'^ ' where d(p is the subdifferential operator of the proper convex lower semicontinuous 

^. , function i^iR'' ^- (— cxd,+oo] and Ct is a second differential operator given by CtVi{x) = 

t"^ ■ ^Tr[a{t,x)a''{t,x)D^Vi{x)]-\- (b{t,x),S/vi{x)), iel^. 

("V^ I We prove the uniqueness of the viscosity solution and then, via a stochastic approach, prove 

OO . the existence of a viscosity solution u: [0,T] x R'' — >■ R'' of the above parabolic variational in- 

>— ' ' equality. 
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a • 1. Introduction 



Viscosity solutions theory was introduced by M.G. Crandall and P.L. Lions [3]. This 
theory allows one to study nonlinear equations which admit as solutions continuous 
functions, without any further smoothness constraints. The classical work in the field of 
viscosity solutions for second order partial differential equations is the survey paper of 
M.G. Crandall, H. Ishii and P.L. Lions [2], where the authors provide several equivalent 
ways to define the notion of viscosity solution, as well as some very general existence 
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and uniqueness theorems. Starting with the study of backward stochastic differential 
equations (introduced by E. Pardoux and S. Peng [10]), generahzed Fcynman-Kac rep- 
resentation formulas have been obtained for the viscosity solutions of semilinear partial 
differential equations. In [4], R.W.R. Darling and E. Pardoux studied elliptic equations 
with Dirichlet boundary conditions; furthermore, in [7], Y. Hu treated elliptic equations 
with homogeneous Neumann boundary conditions. E. Pardoux and S. Zhang extended 
these results in [13] for the case of parabolic systems with nonlinear Neumann boundary 
conditions. 

N. El Karoui et al. [5] considered the case of reflected solutions of one-dimensional 
backward stochastic differential equations, related to an obstacle problem for a parabolic 
partial differential equation. The more general case of backward stochastic differential 
equations involving subdifferential operators and the connection with parabolic varia- 
tional inequalities has been studied by E. Pardoux and A. Ra§canu in [11] and [12]. 

The aim of this paper is to consider the more general case of variational inequalities for 
systems of partial differential equations. The paper is organized as follows. In Section 2, 
we define the extended notion of viscosity solution for a system of parabolic variational 
inequalities. We then formulate the existence and uniqueness result and prove uniqueness. 
In Section 4, we use a stochastic approach in order to prove the existence result; by using 
a backward stochastic variational inequality, we obtain a probabilistic formula for the 
viscosity solution of our system. 



2. Main results 

The goal of this paper is to study the existence and uniqueness of the viscosity solution 
of the following system of parabolic variational inequalities: 

-^{t,x)+Ctu{t,x) + f{t,x,u{t,x))€d^{u{t,x)), te[0,T),xeR^ ^^^ 

u{T,x) = h{x), X6M'', 

where Ctv, with v £ C^(]R'';R'^), is given by 

{CtvUx) ^'^i:Y[a{t,x)a*{t,x)J}^v,{x)] + {h{t,x),m,{x)) 

and T > is the fixed finite horizon. 

We make the following standard assumptions. 

(A.l) The functions 6 : [0, T] x M'' ^ R'' and cr : [0, T] x M'' ^ K''^'' are Lipschitz with 
constant L. 
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(A.2) The function / : [0, T] x M'* x M'^ ^ R'' is continuous and there exists 7 e K such 
that 



.~,l2 



{y - y, f(t, X, y) - /(i, X, y)) <-i\y-y\ 

for ah t 6 [0,r], xGR'', y,yeR''. 
(A. 3) The function h-.M.'^ —^ M.^ is continuous and there exist some Mi > 0, p>0 such 
that, for aUie [0,T], xeR'\ 

\h{x)\<Mi{l + \x\P). 

(A. 4) The function ip-.M.'' — s- (— cxd,+oo] is proper (that is, (p ^ +00), convex, lower 
semicontinuous (Ls.c.) and there exist A/2 > and r > such that 

\(p{h{x))\ < M2{1 + \x\'') \/xeM.''. 

We recall that the subdifferential dip is defined by 

dipiu) = {u* eR'': {u*,v~u)+ip{u) < ^(u),Vw G M'^'}. 

It is a common practice to regard dip as a subset of M.'' x M.'^, by writing {u,u*) g dp 
instead of u* £ dip{u). We define 

Dom(i^) = {ue M'": (p{u) < +00}, 
Dom{dp) = {ue M.'': dp{u) ^ 0}. 

We now give a result which allows us to define the notion of a directional derivative of a 
convex function (see [1] for more details). 

Theorem 1. Let (y3:K'^ — > (— cxd,+oo] be a convex function. Then, for all u G Dom{p) 
and z G M*^ , there exist 



(2) 



, p3{u + tz) - p^{u) ip{u + tz) ~ ipiu) 
LP (u\z):=hTi\ =sup , 

y)(M + tz)-y>(M) . ip(^a + tz)-ip(u) 
ipAu\z) := lim = mt . 

Moreover, the following hold: 

(a) p>'_{u; z) < p'_^{u] z)yu G Dom(v3), Vz G R''; 

(b) p'_{u\ ~z) ^ -p'j^{u]z),\/u G Dom((/3),Vz G M''; 

(c) p'_{u;-) is superlinear and ip'^{u]-) is sublinear; 

(d) if u and z are such that there exists 5 > Q such that u + tz G Dom((^), Vf G {—6, S), 
then (p'_{u;z),p'_^_{u;z) gM. 
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If we take k = 1, then we know that, m every point u £ Doni((^), 

dip{u)=Rn[ip'_{u),ip'^{u)], (3) 

where ip'_ (u) and ip'_^_ {u) are, respectively, the left derivative and right derivative of ip at 
the point u. 

The following proposition generalizes the above characterization to the case of fc > 1. 

Proposition 2. Let ip'.M.'' ^- {—oo,+oo] be a proper convex function and u E Doni{ip) . 
The following statements are equivalent: 

(i) u* G d(p{u); 

(ii) {u*,z)>ip'_{u;z),yzeR''; 
(iii) {u*,z)<ip'_^{u;z)yzeR''. 

Proof. Wc first prove that (i) implies assertion (ii). Let u* G d(p{u). From the definition, 
wc have that 

{u*,y-u)+^{u)<^{y) VyGR^ (4) 

For alH < and any direction z G R*^, we obtain, using (4) with y ^u + tz, that 

{u*,tz) + Lp{u) < ip{u + tz) Vz e M*-', Vt < 

<=> [u ,z)> Vz e M , vt < 0. 

Passing to the limit as t — )• 0, i < 0, wc obtain that 

{u* , z) > Lp'_{u; z) VzgK''. 
Suppose now that (ii) holds, that is, 

{u*,z)>ip'_{u;z) VzGK''. 
From (2), we deduce that 

Consequently 

{u*,tz)<ip{u + tz)-ip{u) VzeM'',V<<0. 

Hence, 

{u*,y-u)<^iy)-ip{u) VyGM^ 

that is, u* € dtf{u). 

The equivalence between (i) and (iii) follows in the same manner. D 
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Let us define, for u e Dom((/3) and z G R'"', 



f*{w,z)— liminf ip_{v;z)^ if '* (u; z) ^ limsup (/3^(u;z). 

veDom{dip} v£Dom{a<p) 

For u S M'^, let (with the usual convention that inf = +00) 

\dip\o{u) ^ ini \dip{u)\. 

If M G Doni{d(p), then there exists a unique u* G K'^, denoted {dif)o{u), such that 
\dip\o{u)^ \{d(p)o{u)\. 

We can now give, using Proposition 2, the definition of a viscosity solution for the 
multidimensional parabolic variational inequality (1). First, we denote the set oi d x d 
symmetric real matrices by S and give the definitions of the superjet and the subjet of 
a function. 

Definition 3. Let u: [0,T] x R'* — 5- R be a continuous function and {t,x) G [0,T] x R^. 
We denote byV'^'~^u{t,x) (the parabolic superjet ofu at {t,x)) the set of triples {p,q,X) G 
R X R'' X 5"^ which are such that for all (s,y) G [0,T] x R'^ in a neighborhood of {t,x), 

u{s,y) < u{t,x) +p{s — t) + {q,y- x) 

+ ^(X{y ~ x),y ~ x) + o{\s - t\ + \y - x\^). 

We similarly define V'^'^uit^x) (the parabolic subjet of u at {t,x)) as the set of triples 
{p,q,X) G R X R'* X S"^ which are such that for all {s,y) G [0,r] x R'' m a neighborhood 
of{t,x), 

u{s,y)> u{t, x) + p{s -t) + {q,y- x) 

+ \{X{y -x),y~ x) +o{\s-t\ + \y - xp). 

Here, r t— > o(r) denotes any function such that linir-j-o -^ = 0. 

Definition 4- Let u : [0, T] x R'' — >■ R''' be a continuous function satisfying u{T, •) — h{-). 
The function u is a viscosity solution of (1) if 

u(t,a;)GDom(.^) V(t,.T) G [0,r) x R'^ 

and if, for all z G R*"'; at any point {t,x) G [0,T) x R'', for any {p,q,X) G V^'^ {u,z){t,x), 
we have that 

p+^Tr{{aa*){t,x)X) + {b{t,x),q) + {f{t,x,u{t,x)),z)>^'^iuit,xy,z). 

Remark 5. li u: [0, T] x R'' ^ R''' is a continuous function which satisfies u{T, •) — h{-), 
then u is a viscosity solution of (1) if and only if 

u{t,x)eDom{ip) V(i,x) G [0,r) xR'' 
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and if, for all 2 G M'', at any point {t,x) £ [0,r) x R'', for any {p,q,X) G 'P'^~{u,z)(t,x), 
we have that 

p+^Tr{{aa*){t,x)X) + {b{t,x),q) + {f{t,x,u{t,x)),z) < >f''*{u{t,x);z). 

We now present the main results. 

Theorem 6 (Existence). Let assumptions (A.1)-(A.4) be satisfied. The multidimen- 
sional parabolic variational inequality (1) then has at least a viscosity solution. 

For the proof of the existence theorem, we use a stochastic approach; wc study a 
backward stochastic variational inequality which allows us to give a probabilistic formula 
for the solutions of (1). We present this approach in the last section. 

Concerning the uniqueness for equation (1), we need to impose more restrictive as- 
sumptions. 

(A. 5) For all u £ Dom(iy9). there exists a neighborhood V oi u such that (9v3)o is 
bounded on Dom((9iy9) n V . 

(A. 6) If M S Dom((p) and 2 G K.'^ such that u + 2 G Dom((^), then there exists a neigh- 
bourhood V oi u such that 

Vv e V^ n Dom(a¥3), 3t > 0: v + tz^ Dom(9v5). 

Theorem 7 (Uniqueness). Let the assumptions (A.l), (A. 2), (A. 5) and (A. 6) be sat- 
isfied. The multidimensional parabolic variational inequality (1) then has at most one 
viscosity solution in the class of continuous functions with polynomial growth. 

3. Proof of the uniqueness theorem 

We first prove the following lemma. 

Lemma 8. If ip:M.'' — > (— cx),+oo] is a convex function, proper and l.s.c., then we have: 

(a) for every u,v G Doia{ip), 

ip'_{u,u — v) >ip'_^{v,u ~ v); (5) 

(b) for every w,w G Dom{d(f) and 2 G M*^, (p'^{u;z) < (p'_{u;z) and (p''*{u;z) > 
f'+{u;z); 

(c) for every z € M.'' , (p'^{-;z) is l.s.c. on Dom((/3) and if''*{-; z) is upper semicontinous 
(u.s.c.) on Dom((^); 

(d) for every u,d G Dom((/3), if 

limsup (inf [|a(/3|o(w' - t{u - v)) + \dip\o{v' + t{u - v))]] < +00, (6) 



264 Maticiuc, Pardoux, Rd§canu and Zalinescu 

then 

Lp'^{u,U — v) > ip''*{v,U — v). 

Proof, (a) By the convexity oi tp, for t G (0, 1), we have 

i^((l - t)u + tv) + ip{tu + (1 - t)v) < ip{u) + lf{v). 

Therefore, 

ip{u — t{u — v)) — ip{u) ip{v + t{u — v)) — ip{v) 

t - t 

and passing to the hmit as t \j 0, wc obtain ip'_ {u, u — v)> (p'^{v, u ~ v). 

Results (b) and (c) follow in a standard way from the definitions of ip'^ and ip''* , so we 
skip their proofs. 

(d) For simplicity, let us define z = u — v. From (6), we can find M > and r > such 
that, for u' E B{u, r) D Dom{dip), v' € B{v, r) D Doni{dip), there is a t^'y < such that, 
for t G {tu',v',0), we have 

Idifiloiu' + tz) + \dLp\a{v' - tz) < M. 

Let m', v' and t be as above. Then, 

if{u' + tz) - ip{u') _ Lp{u' + t{u' - V')) - Lp{u') ip{u' + tz) - ip{u' + t{u' - v')) 

i ~ i t 

ip(u' +t(u' — V')) — (p(u') ,,„,,, s , , ,^^ 

> VA ^ ^ )) ^\ J ^ ((5<^)q(^/ + i^)^ z-{u'~ v')). 

It follows that 

^{u' + tz) - ip{u') ^ ipju' + t{u' - v')) - ^{u') 

> Al\z — {u —v)\. 

t - t I y Jl 

Passing to the limit as t /^O, we have 

ip'_{u'; z) > ip'_{u';u' - v') - M\z - {u' - v')\. 
In a similar manner, we obtain that 

f'+iv'; u' - v') > ip'+{v'; z) - M\z - {u' -v')\. 
Combining the last two inequalities with (5), we deduce that 

ip'_{u']z)>Lp'^{v']z)-2M\z-{u' -v')\. 
Passing to the limit as m' — ?► m and v' — ?> w, we obtain that (^K"! z) ^ ^''*(v; z). D 
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Remark 9. From point (c) of the above lemma, we can consider in Definition 4, 
'P'+{u,z) (resp., :^^-(i;,z)) instead of 7'^'+(u,z) (resp., 'P'^~{v,z)), where ■^'^^+ and 
P^'~ denote the respective closures of operators p^'+ and V^'~ (see the definitions in [2]). 

Proof. Let u,v € C([0,r] x M'^;]R'^) be viscosity solutions for (1). Wc must prove that 
u = w on [0,T] xR''. 

For some A > to be made precise below, wc carry out the following transformations: 

u{t, x) = c^*C^{x)u{t, x), v{t, x) = e^*^-i (a;)u(t, x), 

where ^(x) = (1 + \x\'^)'^^^ , with /i larger than the order of growth of u and v. 
Then, u and v are bounded and they are solutions, in the viscosity sense, of 



' du{t,x) 



+ Au{t, x) + /(i, X, u{t, x)) e e^*C"i (.t) dip{e-^^£,{x)u{t, x)), 



(7) 



dt 

te[0,r),.TeR'*, 

Vu{T,x)^e^^i-^{x)h{x), cceR^ 

where, for i = 1, fc, 

{:^tu)^{t,x) = {Ctu),{t,x) + ((aa*)(t,a;)r'(a^)Ve(x), Vu,(t,a;)) 

and 

fit, X, u) = e^*r' (2;)/(i, •^, e-^*C(-^)u) " AS 

+ iTr[(aa*)(t,a:)r'(x)D2e(a;)]fZ+(6(i,.T),r'(a:)VC(a;))w. 

Since lim|2.|_j.^QQu(i,x) = lim|3;[_j,+ooS(^7 2;) = and u{T,x) = v{T,x), there exists 
{to,xo) G [0,r) X R'' such that 

e:^\u{to,xo)~v{tQ,XQ)\ > \u{t,x)-v{t,x)\ V(t, x) G [0,r] x R'^. 

Let us suppose that u^v. This implies that ^ > 0. We set 

zo = T(u(^o,a;o) - v{to,xo)). 
Then, 

0:^ {u{ta,xo)~v{to,xo),zo) > {u{t,x) ~ v{t,x),zo) V(i,2;) e [0,T] x R'^. 
We let, for a > 0, 

<^a{t,x,s,y) :=(u(t,x) -v{s,y),zo) - -{\t - s\'^ + \x ~ y\'^). 
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Since lunswpi^i^i i^^^^ait,x,s,y) < and $0(^0,2:0:^0,2^0) = ^ > 0, there exists 

{ta,Xa,Sa,ya) G ([0,T] X K'^)^ SUCh that 

Ma ■■=^a{ta,Xa,Sa,ya) =SUp$Q. 

By Lcnniia 3.1 from [2], we have: 
(a) 

Km a{\xa - Val'^ + \ta- Sa\'^) = 0; 

(b) whenever (i,x) is an accumulation point for {ta,Xa) as a — > +00, wc have that 
hm Ma = {u{i, x) - v{i, x),zo) = 0. (8) 

Since, for large a, {ta,Xa) remains in a compact subset of [0,r] x 'R'^, there is at least 
one accumulation point {t,x) for {ta,Xa) as a — > +00. We can suppose, without loss of 
generality, that (ta,Xa) — >■ it,x). Of course, from (8), t <T, and we can also assume 
that ta,Sa G [0, T) for every a. Another consequence of (8) is that u(t, x) — u(i, a;) = Ozq. 
Indeed, 

\u(t,x) — v{t,x)\ <0 = {u{t,x) — v{t,x),zo) < \u{t,x) — v{t,x)\, 

from which we conclude that u{t,x) — v{t^x) = |'«(i,a:) —v{t, x)\zq = Ozq. Let us now apply 
now Theorem 3.2 from [2], which asserts that, for every a > 0, there exist X,Y € S*^ such 
that: 

(a) {a{ta - Sa)-,a{Xa -ya),Xa) E P'^'^ {u, Zn){ta, Xa)] 

(b) {a{ta - Sa),a{Xa -ya),Ya) G ^^'^ (w, ^o) (Sa , J/c) ; 

From the definition of the viscosity solution and Remark 5, we have: 

>e^*°r^(2:a)¥'Uc-^*°e(2:a)S(ia,2;„);zo), 

a{ta -So) + £'a(x„-y^),Y^{Sa,ya) + {f{Sa,ya,v{Sa,ya)),Zo) 
<e^'-C\ya)v''*ie-^'''(iya)viSa,ya)]ZQ), 

where, for qeR'^, X eS'^, 

4x(t, 2:) = i Tr[(aa*)(t, x)X] + {bit, x), q) + ((aa*)(t, x)^' {x)Vax), q). 
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Subtracting the two inequalities, we obtain 



>e^*°^-^(xa)(pt(c-^*°^(a;„)u(i„,x„);zo) 

From point (c) of Lemma 8, we get 

\immie^*''£,~^{xa)'p'Ae~^^''£.{xa)u{ta,Xa);zo) 

limsupe^''°^"^(ya).y9''*(c"^'*°C(ya)w(sa,ya);2o) 

a— J-+00 

<c^'rHxW'*{e-^'^ix)viixy,zo). 

Now, by point (d) of Lemma 8 and assumptions (A. 5) and (A. 6), 

^:(c-^*^(£)fi(£,x);zo) > v''*{e-^'axMi,x);zo) 
since u{t,x) — v(t,x) = 9zq. It follows that 

lirRmi[c^*-C^{Xc,)ip'^{c-^*^^{Xa)u{ta,Xa);Zo) 

a— ^+oo 

-e^^°r'(2/a)¥'''*(e-^^°C(ya)«(Sa,2/a);2;o)] 

>e^'r'{iWAe-^'aiMi,xy,zo) 

-e^'r\S:W'*ic'^'ax)v{i,xyzo)>0. 
On the other hand, we have that 

= ^Tr[aa* {ta, Xa)Xa -~a-a*{sa,ya)Ya] 

+ a{b{ta,Xa) -b{Sa,ya) 

+ aa*{ta,Xa)^~^{Xa)y^iXa) - tTCT* (Sa , J/q)^"^ (?;«) V^(ya), a;^ - J/a) 

< ^a\a{ta,Xa) - a{sa,ya)\'^ 

+ a\Xa -ya\i\b{ta,Xa) ~b{Sa,ya)\ 

+ \(T(T*{ta,Xa)C^{Xa)'^^iXa) - Cr<J*{Sa,ya)C'^{ya)'^^{ya)\) 

< 3a{Ll + Li){\Xa - 2/qP + \ta - Sap), 



(9) 
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where Li is the Lipschitz constant of (cr, 6,CT(7*^^^V^). Hence, 

a— > + oo 

Finally, from assumption (A. 2), 

lim {f{ta,Xa,u{ta,Xa)) - /(Sq, I/q , i^(Sa , 2/q)), Zq) 

a— f+oo 

= (/(t, X, u{i, x)) ~ f{i, X, v{i, x)),zq) 
^{-X+^TT[aa*{i,x)r\i)^'ax)] + {b{i,i),r\x)\/ai))W 

+ e^'r\x) (/(£, X, c-^'^xMi, x)) - fii X, c-^'ax)vii x)),zo) 
< {-X+^Tiiaa*{i,x)r\x)^'ai)] + (6(t,i),r'(i)VC(i)) +7)^- 
By taking 

A>7 + sup[iTr[aa*r'D'e] + (&,r'V0], 
we get 

lim {f{ta,Xa,u{ta,Xa)) " /(Sq , 2/q , w(Sq , J/a)) , ^o) < 0. 

Passing to the limit as a — >■ +00 in (9), we obtain a contradiction. Hence, u = v and 
so u = v. 

D 

4. Proof of the existence theorem 

In this section, we prove Theorem 6. This is obtained via a stochastic approach. Using a 
certain backward stochastic variational inequality, we will obtain a generalized Feynman- 
Kac representation formula for the viscosity solution of (1). 

Let {Wt'. i > 0} be a d-dimensional standard Brownian motion defined on some com- 
plete probability space {il,J-,¥). We denote by {J-j: t>0} the natural filtration gener- 
ated by {VFt: t>0} and augmented by J\f, the set of P-null events of T: 

T^ = a{Wr: 0<r<t}VA/'. 

4.1. Backward stochastic variational inequality 

We consider the following backward stochastic variational inequality 



dY,+F{s,Y,,Zs)dsedip{Ys)ds + Z,dWs, se[o,r]. 



(10) 



where: 
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(I) F-.nx [0, T] 



pkxd V itj) fc 



satisfies, for some a € M, /5 > 0, 



(i 
(ii 
(iii 
(iv 

(v 



F{-,-,y,z) is a progressively measurable stochastic process (p.m.s.p.), 

y I — > F{uj, t, y, z) is continuous, a.s., 

{y - y',F{t, y, z) - F{t, y', z)) < a\y - y'\^, a.s., 

\Fit,y,z)-F{t,y,z')\<L\\z-z'lei.s., 



e( I F*{t)dtj <oo, 

Vi?>0, where i^*(i):=sup{|i^(i,y,0)|: |y| < i?}. 



(11) 



for aU te[0,T], y,y' eR'', z,z'eM'=^'*; 
(II) if.SJ' — > (— 00,+cxd] is proper, convex and l.s.c; 
(III) the terminal date ^ G L'jr^{rt;'R'') is such that 

E(^(e)) < ^. 

For the proof of the next theorem see [9] . 



(12) 



Theorem 10. Let the assumptions (I)-(III) be satisfied. There then exists a unique triple 
{{Yt,Zt,Ut): te[0,T]} of p.m.s.p. such that: 



sup |y,p+ f {\Zs\^ + \Us\^)ds 

s£[0,T] Jo 



(a) E 



(b) E/ ip{Y,) ds < +oo; 
Jo 

(e) {Yt,Ut)€dip, 



< +oo; 



)dt a.e. on fix [0,r]; 



(d) yt+ f Usds = ^+ f F{s,Ys,Zs)ds- f ZsdW^ for all t>0, a.s. 
Jt Jt Jt 



Let 



J"* = a{Wr ~Wt: t<r<s}yj\f. 



From assumption (A.l) of Section 2, it follows (see, for example, [6, 8]) that for each 
{t,x) £ [0,r] X M'* there exists a unique continuous {J-"*} -p.m.s.p. {Xl'^: s G [0,T]} 
solving the stochastic differential equation 



Xl'^ = x+ f b{r,Xl:'')dr+ j a{r,Xl:'')dWr. 
Jt Jt 



(13) 
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Remark. Note that, in this section, one can relax the assumption (A.l) by requiring 
only the continuity of b and cr in t € [0,T], instead of Lipschitz continuity. 

The following proposition summarizes some well-known properties of the solutions of 
equation (13). 

Proposition 11. Under assumption (A.l), the solution {Xl'^: s G [0,r]} of equa- 
tion (13) satisfies that for all p>2, there exists some constant C = C'{p,T,L) > such 
that, for all t,i £ [0,T], x,x e R'^ : 

(i) E sup \Xl'-\P<C{l + \xn 

se[o,T] 

(ii) E sup \Xl'''-Xl'^\P<Cil + \x\P + \i\P)i\x-S:\P + \t-i\P/^). 
se[o,T] 

From Theorem 10 with ^ = h{Xj,^), F{u!,s,y,z) ~ f{s,Xl'^{uj),y), it follows, under 
assumptions (A.1)-(A.4), that for {t,x) G [0,r] x M'', there exists a unique triple of 
{J^]}-p.m.s.p., {(Fg'^, Zl'^', t/*'^): s £ [0,T]}, solving the backward stochastic variational 
inequality 



y;^"+/ C/;-"dr = /i(X^")+ / /(r,x;■^F;^")dr 

for ah s G [t, T] a.s. 



Zl.'"" dWr 



(14) 



t,X ryt,: 



with (r;-^, t/*'^) e av', P X As a.c. on rj x [i, T], and Yy-' = F/'"^, Z\ 
Moreover, we have the following properties of the solution of (14). 



[/ 



:OVse[0,i]. 



Proposition 12. Under the assumptions (A.1)-(A.4), we have that there exists some 
constant C > such that for all t,i G [0,r], x,x G W^: 

(i) E sup |y,*'n^<C(l + |a;|2); 
se[o,T] 

(ii) E sup |y^*.--yp|2 

se[o,T] 



<E 



\h{X*^^) - h{X*/)\^ 



+ / |l[t,T](r)/(r,X*^^r/^-)-l[,;^](r)/(r,X*■^y;■-)|^dr 



We define 

u{t,x) =¥*''', {t,x)G[0,T]xR'^, 

which is a determinist quantity since Y^ '^ is J-^ = cr {A/"} -measurable. 



(15) 
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From the Markov property, we have that 

u{s,Xr)=Y:^^. (16) 

For the proof of the next proposition, see [11]. 

Proposition 13. Under the assumptions (A.1)-(A.4), the function u defined by (15) 
satisfies: 

(i) u{t,x) e Dom((^),V(t,.T) G [0,r] x W^; 

(ii) |w(t,a;)|<C(l + |.T|P),V(t,a-)e[0,r]xR'^; (17) 

(iii) ueC([0,r] xR'';R'=). 

4.2. Connection with variational inequalities for systems of 
partial differential equations 

We prove that the function u defined above is a viscosity solution, in the sense of Defi- 
nition 4, for the muhidimensional parabohc variational inequality (1). 
Theorem 6 is a consequence of the following result. 

Theorem 14. Let the assumptions (A.1)-(A.4) be satisfied. The function u defined 
by (15) is then a viscosity solution for the multidimensional parabolic variational in- 
equality (1). 

Proof. Let z ^M.^ , {t,x) G [0,7] x R'' and {p,q,X) e 'P'^'+{u,z){t,x). Let us define, for 
simplicity, 

V{t,x,u,p,q,X)^p+^Tii{aa*){t,x)X) + {b{t,x),q) + {f{t,x,u),z). 

Suppose, contrary to our claim, that 

V{t, X, u{t, x),p, q, X) < ifi'^ {u{t, x);z). 

It follows, by the continuity of /, u, b, a, that there exist e > 0, S > such that for all 

\s-t\<d,\y-x\<S, 

V{s,y,u{s,y),p + e,q+{X + eI){y-x),X + eI)<'f',{u{s,y);z). (18) 

Now, since {p,q,X) G V'^'^{u,z){t,x), there exists <6' <d such that 

{u{s,y),z)<'i{s,y) (19) 

for ah se[0,T], s^t,yeW\y^x with 0<s-t<S', \y - x\ < 5', where 

*(s, y) = {u{t, x),z) + {p + e){s - t) + {q,y - x) + ^{{X + el){y -x),y-x). 
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T = mi{s>t: |X*'^~x|>5'}. 
^Y^t,.jt,.^^^^Y:'^^z),{Zr,z)), t<s, 



Y^^- = {u{T,Xl'%z)+ j ((/(r,X;^%w(r,X;'-)),z)-([/;^%2))dr-y Zl^^ AWr 

{Ul--,z) 6 [^'_(r,*.-z),^V(y,*'-z)]. 

From Ito's formula, it follows that 

(i;*'%Z^-):=(vI/(s,X*^-),(Vvl'a)(,s,X*^-)), t<s<t + 6\ 

satisfies 

"9* 



!;*'" = *(t,x*'=^)- f 



dr 



(r,x;--) + /:,*(r,x;- 



dr 



ZT dWr 



Let (y;*■^ zl'-) = (n*-- - n*>-, z*'- - z^-). 

We have 






(20) 



dr 



Zb'^dW^. 



We note that, from Lemma 8, point (c), 

{Ul'\z)>^'_{<s.Xln;^)>^:{n{s,Xl'-)-z), P^dia.e. 
Moreover, the choice of 6' and r implies that 

vl/(T,X^-)>Kr,X*-),z). 
From (18), it follows that 

^',{u{s,Xl^-)-z)>V{s,Xl'\u{s,Xl'-),p + e,q + {X + eI){X'f~x),X + eI) 



'd7 



{s, X^) + £,*(«, X^) + {f{s, Xr, u{s, X^)), z). 



These inequalities and equation (20) imply that y^ '^ > or, equivalently, 

*(i,a;)>(u(i,a;),z). 
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which contradicts the definition of ^. Hence, we have 

V(t,x,u{t,x),p,q,X) > ip'^(u(t,x);z). 
This proves that u is a viscosity solution of (1). D 
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